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A two-level system is considered which may perform flip-processes by coupling to a classical bath 
represented by a stochastic field. The time evolution of the density matrix leads to a stochastic 
equation with a multiplicative noise. The corresponding Fokker-Planck-equation (FPE) for the 
probability density depends on the matrix elements of the underlying density operator. The solution 
of the FPE is parametrized in terms of a new conserved quantity a, which follows from the stochastic 
equation. The parameter a is interpreted as a measure for the persistence of quantum information. 
The FPE exhibits an single unique stationary solution which is different from Boltzmann's law. The 
eigenvalues of the time dependent part are calculated exactly, leading to discrete relaxation times 
characterized by two quantum numbers and the ratio of Planck's constant and the coupling strength 
to the bath. The entropy S{a) is analyzed as function of the quantum number a. In case of a = 1 
the system is in a pure state whereas for a ^ 1 a. mixed state is realized. If a second two-level 
system is included, immersed in the common bath, both noninteracting two-level systems become 
mutually entangled. The annealed entropy is in that context non extensive. 

PACS numbers: 03.65.Yz, 03.67.-a, 05.70.-a, 03.65.Ca, 03.65.Ud 

There is an increasing interest in the role of macroscopic environments to our understanding of the basics of quantum 
theory. The knowledge of the implications of the quantum theory to other theories, especially to the statistical 
mechanics and the domain of validity has captivated scientists from the beginning of quantum description. In such a 
context, the presence of an environment is commonly thought as entanglement, decohering and mixing properties of 
quantum system. Generically, an environment is assumed to be a noisy reservoir or a heat bath. Whereas in common 
interpretation of statistical mechanics the heat bath is unspecified, in quantum systems a heat bath can also provide 
an indirect interaction betv^^een otherwise totally decoupled subsystems and consequently a means to entangle them 
Q, 1^ . In simple example for the entanglement between two qubits due to the interaction with a common heat bath 
has been explicitly shown in 0|. Whereas in that paper the bath is described by a collection of harmonic oscillators, 
it seems to be more reasonable to specify the bath by stochastic forces represented by stochastic fields. From a more 
general point of view we expect the bath should be better described in a stochastic manner and not by deterministic 
forces. In the present paper we consider a two level system (qubits) which are able to perform fiip processes by a 
coupling to classical stochastic fields. Thus we bridge the gap between quantum and classical probability theory. This 
problem is related to many other questions of quantum optics and quantum electronics where quantum statistical 
aspects arising from the intrinsic quantum character of the system while the possible time-dependence of system 
parameters may be interpreted as the influence of classical thermal fluctuations. 

A very simple idealization of many quantum systems is a single two-level system which can exist in either of two 
states, |-f 1) and | — 1). The quantum mechanics can be represented by an U{2) matrix representation with the mapping 
— 1) = (0,1) and 1+1) = (1,0) and the Pauli operators (T^,cr+,cr^ in their standard representation. Obviously, cr+ 
lifts the two level system from the ground state | — 1) to the excited state |-|-1) while a~ drops the system to the 
ground state. The Hamiltonian of the two level system in an external time-dependent complex scalar field E{t) can 
be written as 

H = i - ^hr] [cr+E{t) + (J-E*{t)] 

where the energy difference between the ground state and the excited state is given by hcu and the coupling between 
the two level systems and the external field is defined by rj. The Hamiltonian allows the explicit formulation of the 
time evolution of the density operator g, which is formally given by the von Neumann's equation. The diagonal matrix 
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elements satisfy the following equation 



'^%r " ljivE{t)Q-+-l^vE*{t)g+., (1) 
whereas the non-diagonal element obeys 

z^i— ^ = fiujg+- - -h-qEit) [g — - g++] . (2) 

The other element follows from g^ = The four components of the density matrix are not completely indepen- 
dent. The sum of the diagonal elements + g = const corresponding to the necessary normalization condition 

trg = g+-i- + = 1. The external field E{t) represents the influence of a classical environment on the quantum 
mechanical system. We assume, that this field may be expressed by a complex white noise, E(t) = Eg [£,i{t) + 1^2 (i)], 
with ^1/2 (i) corresponding to normalized Wiener processes £,i/2{t)dt = dWi/2{t) with the noise strength Eq. Thus, the 
density matrix g describing the probability distribution of the quantum process is itself a randomly varying quantity. 
On the other hand, the system of evolution Eqs. I|1I2|) is an ordinary set of Langevin equations. In other words, it 
should be possible to map the dynamics of these equation onto a Fokker-Planck equation. As argued above we discuss 
the dynamics in terms of the Fokker-Planck equation (FEE) . To this aim, the Langevin-equations should be reduced 
to a set of independent quantities representing the four components of the density matrix. The first step is suggested 

by the normalization condition and leads to the parametrization = cos^ ^ , g = sin^ ^ , g^ = ii?e*"^ . 

Therefore the set of the three independent Eqs. |Q][2J) is parametrized by the angles < < tt and < $ < 27r as well 
as the real function R{t). In terms of these new variables Eqs. 1^121 can be rewritten as three independent stochastic 
equations of the form 

sinOe = 77£:oi? [cos 0^2(0 -sin 0^1 (<)] 

R — T] Eq [cos (1)^2 {t) — sin (j)£^i{t)] cos <d 
R<j) ^ ~ujR ~ rjEo cos 9 [sin 0^2 (t) + cos <p^i (t)] . (3) 

From the first and second relation follows immediately the conservation law 

cos^ e{t) + R^{t) ^c? ^ const. (4) 

This conservation law is valid for each trajectory. In the following studies we choose a = 1 which corresponds to a 
pure state which will be discussed in detail below. We may use Eq. as a constraint in order to eliminate the second 
relation in Eq. ©. For a = 1 Eq. Q suggest to set R(t) — sin9(i). It results a coupled system of two stochastic 
differential equations 



de = ^EQ [cos<\)dW2(t) - s\n<\)dWx(t)\ 



d<t) = -ujdt - [sin 4>dW2{t) + cos (t>dWi (t)] , 

tan B 

which may be interpreted in the sense of the Ito calculus. From here, we get by application of the standard formalism 
[3 the FPE 

^F(e,0,<)-£p(e,(/.,i), (5) 

where the operator L is defined by 

" 902 + 4tan2e902 + ' 

The noise strength E^ is incorporated in the coupling parameter 77. The last equation describes the time evolution 
of the probability density P(B, (j)) to find a quantum statistical density operator of a stochastically driven two level 
system. The density operator is parametrized by the two remaining degrees of freedom and and the initial 
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condition is given by P(6, 0, to) = 5 {Q — Go) S {(j) — (j>o) , where the initial values 9o = 6(to) and (j)o = <p{to) may be 
directly obtained from the initial quantum mechanical state. Without any coupling between the quantum mechanical 
two level system, ij — 0, the FPE is reduced to the simple Liouville equation 

with the well known solution 

P(e, ^,t) = 5{e- Go) (5 (0 - 00 + |i) 
This result corresponds to the dynamics of the density matrix 

Q++{to) g+^{tQ)exp{-iujt/2}\ 
[g-+{to)exp{iu;t/2} e__(io) J 

with fixed diagonal elements but periodically time-dependent non-diagonal components. In order to solve the general 
Fokker-Planck equation we introduce the new variable (p = (p + ujt/2. This leads to the elimination of the drift 
term in Eq. (jSJ 

dt ' 4tan2e9(^2 ^ > 

with P(8, (f>, t) = P(6, (j) + t). The solution of these equation can be done by an expansion in terms of eigen- 
functions. To this aim we use the ansatz 

P(e,(p,t)=Q(e,(p)exp{-AO 

and obtain the eigen-value equation 

92 1 92 \ 

with A ~ ifK/A:. The parametrization of the density matrix requires the periodic behavior Q(0, (p) — Q(0, </? + 27r). 
Therefore, the eigen functions Q(0, (f) can be written as 

0(9, (^) = 5(e)exp{im(^} 

The transformations 5'(8) = y (cos Q) sin^/2 q a^cj subsequent substitution x — cos Q leads to the equation 

A - 1 1 + m2 



2\ ,1 n , , /m2+A-l l + m2 \ 



The two independent solutions of the last equation arc given in terms of the associated Legendre functions P^' (x) and 
Q'i{x). It results 



The eigenvalues A follow from the constraint, that the probability density P(8, (/)) should not offer a singular behavior. 
Here, this means (1 — a;2)i/4 yf^^^-j should not singular for x ±1, which leads to 



A„„ = {2n - i f + 2{2n - 1)^1 + m'^ + 1 V?i e N. (8) 

The knowledge of the eigenfunctions and corresponding eigenvalues enables us to construct the initial distribution and 
consequently, to find out the time evolution of the system. Since all eigenvalues A„„j has to be positive semi-definite, 
the asymptotic distribution Poo(0,(/3) = P{<d,(p,t — > 00) belongs to the eigenvalue A = 0. It is easy to check, see 
Eq. ijSJ that this is relalized by setting n = m = 0. There exist only one eigenfunction with Aqo = 0, namely Poo (6, 4>) 
= const, for all Q and (p. On the other hand, the final distribution function Poo(6,<?ii) defines a stationary state. 
Thus, Poo (6, (/)) gives rise to a constant probability current, namely 

J _ 2 9Poc J _ dPpo ^ uj ^ 

® ""^ 99 ' ^" 4tan2e 86 ~^ 2 °" ' 
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These equations has the solution 



2 

Poo(0, <j>) ^ —J4, = const. 



2 I i?e-'* 1 - cose ) ■ 



and the constraint Je = 0. The probabihty distribution function P(0, (f>, t) describes the statistical properties of the 
quantum system. The variables of the function P are the matrix elements of the underlying density operator of the 
quantum model. In order to interpret this on a first view curious result, let us go back to Eq. (0J. Especially we are 
interested in the meaning of the new quantum number a in more detail. In the general case the density matrix reads 

1 / 1 + cos e i?e** 

2 I 

So it results det p = ( 1 — ) /4 and trp^ = (l + a^)/2. From here we conclude that a is a real quantity with < a < 1 . 
Obviously, the before analyzed case a = 1 corresponds to a pure state, tr^^ = 1 and = g, respectively. Apparently, 
a 1 defines a mixed state. Remarkable is that the degree of mixing, expressed by a 7^ 1, offers a constant value 
independently from the time evolution of the external field E(t). Insofar a quantum state reveals a kind of persistence 
under the influence of a classical bath. The curiosity of the result consists of that the probability distribution of the 
quantum object is again subjected to a statistics, i. e. the quantum probability is itself an object of a statistics. 
Regarding this situation the preparation of the quantum object will be essential, since the information, included in 
the quantum statistical operator of a closed quantum system, is determined exclusively by the state of the system 
after its preparation. Notice that (p and can be geometrically interpreted as spherical coordinate. In case of an 
arbitrary < a < 1 in Eq. the angle Q is restricted to a region Qq < Q < tt — Qo with cos0o = a, see Eq.Ql. 
This restricts to a strip around the equatorial cycle of a sphere. 

Let us discuss the entropy as a measure of information about the two level system. Using the density matrix obtained 
in Eq. Q and its representation in the form g = al + 2b- a with appropriate chosen scalar a and vector b as well as a 
representation of an arbitrary function /(al + 2b- a) in terms of the Pauli matrices a ^2Qj. we get the entropy, defined 
as 5 = — tr^ln^: 

5* = i [21n2- (1 + a)ln(l + a) - (1 - a) ln(l - a)] . 

That means the complete quantum information of the two level system is independent from the time evolution of 
the classical field which drives the system. The background for this apparently surprising result is the understanding 
of the measurement processes applied to the determination of the density matrix g. In the present case the system 
follows its evolution equation and the entropy is completely defined by the initial preparation by which the degree of 
mixing a is fixed. A further determination of g at a later time is not necessary. Thus, the density matrix can here 
interpreted in the sense of Bayesian statistics as a certain degree of believe suggested by the initial preparation. The 
traditional physical interpretation of a probability distribution function is the Gaussian frequency concept which is 
reflected in the ensemble theory. Hereby, the density matrix is measured again at ig- Since a probability function bases 
on a frequency distribution, the determination of g{te) requires multiple measurements with identical experiments. 
As a consequence, we get now an averaged density matrix {g{te)) which includes also the average procedure over the 
external field, i.e. 

idite)) = I ?(e, 0)P (e, 0, ie) dQd^ 

and the actual informational content is now given by 

Sa = ~tY{g{te)\n{g{t,)), (10) 

which is equivalent to the thermodynamic definition of the entropy in case of an equilibrium state. The index a stands 
for annealed, so below. With other word, the classical and quantum average procedure do not commute provided the 
the mixed state of the initially realized state can be characterized on the basis of quantum mechanics. 
From this point of view, we have to distinguish between the quenched average with respect to the external stochastic 
field and the annealed average. The last case is realized by the average procedure made in Eq. l{Tn|l while the quenched 
case is defined as Sq = — (tr^ln^). In that case we obtain 

Sg = -1 [21n2- (l + a)ln(l + a) - (1 - a) ln(l - a)] , 
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where the average is reahzed by means P{Q, 4>, t). The difference between both expressions Sa and 5*^ is a consequence 
of the assumed non-equihbrium character of prepared initial state. If we start directly from a system with a maximum 
mixing state, a = 0, and therefore from the equilibrium state, we get always Sq — Sa- The discriminate between the 
quenched and annealed average procedure is not necessary, if we analyze quantities which are linearly dependent on 
the density matrix, for instance the averaged spin S^q — trgS^ = trgS^ — S^a ■ Thus, the measurement of primarily 
non-thermodynamical quantities like internal energy or averages of spin operators is not affected by the commutation 
of the classical and quantum mechanical average, while all non-linear functionals of the density matrix, like degree of 
mixing or informational content depend strongly on the order the applied average procedures. Now let us generalize 
the model to two spins in the same classical field. Following the procedure as before the FPE for P(6i, 62, 4>2t) 
has the same form as in Eq. ©. However where the evolution operator reads L = Li + L2 + L12. The additional 
term Li2(6i, 62, (^2, t) is originated by the presence of a common bath, represented by the same stochastic field, 
for both spins. As a consequence we find 

(^(61, 0i) 0^(62,02) {d{Qi,<f>i)) ® {d{Q2,^2)) . 
Whereas the total quantum entropy S'tot = 81+82 is additive, the total annealed entropy, defined as 



Stot.a = tre(ei,0i) (8)^(62, 02) In e(ei,0i) «)?(62,(?!'2 
is non extensive. A short time expansion in an interval At leads to 



m = p{to) 



71 At 

tanBi tanB2 



Although the two spins are completely without mutual interaction the system reveals a correlation between both 
particles. The quantum spins become entangled due to the common classical bath. A complete incoherent quantum 
state becomes coherent due to the coupling at a common bath. Concluding the paper we have considered a quantum 
system immersed in a classical bath represented by a stochastic field. In this sense the matrix elements of the density 
operator are subjected to a statistical description. They are the variables in an underlying FPE. The solution of the 
FPE is parametrized in an additional conserved quantity a by which mixed and pure states can be distinguished. 
When a second two-level system is immersed in the same classical bath, the initially decoupled quantum systems 
become mutually entangled. Insofar correlated baths violate the rules of the classical concepts of thermodynamics. 
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